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\ 1 INTRODUCTION 



ABSTRACT 

We report the first results of our programme to obtain multi-epoch radial velocity 
measurements of stars with a strong far-UV excess to identify post common-envelope 
binaries (PCEBs). The targets have been identified using optical photometry from 
SDSS DR4, ultraviolet photometry from GALEX GR2 and proper motion infor- 
mation from SDSS DR5. We have obtained spectra at two or more epochs for 36 
targets. Three of our targets show large radial velocity shifts (> 50kms _1 ) on a 
timescale of hours or days and are almost certainly PCEBs. For one of these targets 
(SDSS J104234. 77+644205. 4) we have obtained further spectroscopy to confirm that 
this is a PCEB with an orbital period of 4.74 h and semi-amplitude K = 165kms~ . 
Two targets are rapidly rotating K-dwarfs which appear to show small radial velocity 
shifts and have strong Call H+K emission lines. These may be wind- induced rapidly 
rotating (WIRRing) stars. These results show that we can use GALEX and SDSS pho- 
tometry to identify PCEBs that cannot be identified using SDSS photometry alone, 
and to identify new WIRRing stars. A more comprehensive survey of stars identified 
using the methods developed in this paper will lead to a much improved understanding 
of common envelope evolution. 
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Post common-envelope binary stars (PCEBs) can be loosely 
denned as short period binary stars containing a compact 
object. Cataclysmic variable stars (CVs), in which a white 
dwarf accretes matter from a low-mass star that fills it Roche 
lobe, are one example of PCEBs. Compact objects are gen- 
erally the result of stellar evolution involving at least one 
phase when the star is a red giant, but the size of a giant 
star exceeds the current separat ion of the s tars in a PCEB 
by a few orders of magnitude. iPaczvnskH (|l976l ) outlined 
the solution to this puzzle in a short paper on the origins of 
CVs that describes the main features of what is commonly 
known as common envelope evolution. In the case of CVs, 
the evolution begins with a star of a few solar masses with 
a less massive companion in a wide orbit. The evolution of 
the more massive star through the red giant phase is ter- 
minated prematurely when it exceeds the size of its Roche 
lobe. The resulting mass transfer onto the low-mass com- 



panion is dynamically unstable, so the companion is unable 
to accrete the material from the red giant. The material 
forms a common envelope around the dense core of the red 
giant and the low-mass companion. If the companion is not 
massive enough to force the envelope to co-rotate then the 
evolution proceeds on a dynamical timescale as dynamical 
friction forces the low-mass companion to spiral in through 
the outer layers of the red giant. The orbital energy lost by 
the companion is transfered to the outer layers of the red 
giant, which are ejected to reveal the hot, degenerate core 
of the red giant and, if it survives, the low-mass companion. 
Angular momentum loss by a combination of a magnetic 
stellar wind and gravitational wave radiation, or the evolu- 
tion of the companion may lead to a second phase of Roche 
lobe overflow, this time from the low-mass companion onto 
the white dwarf, i.e., the formation of a CV. For this rea- 
son, detached binary stars containing a low mass star and 
white dwarf with a n orbital period of about a da y or less are 
known as pre- CVs (|Schreiber fc Gansickdliool h 
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Figure 1. The GALEX FUV and NUV response (arbitrary units) 
compared to IUE spectra of 3 white dwarfs (HD 149499 B DOZ1; 
G191-B2B, DAI and EG 274, DA2), and an active solar-type star 
EK Dra (in order of decreasing FUV flux). All fluxes scaled to a 
distance of 10 pc. 



The compact object in a PCEB may be a neutron star 
(NS), black hole (BH) or white dwarf (WD). Short period bi- 
nary stars containing hot subdwarf stars (sdO or sdB stars) 
are also considered to be PCEBs. The companion star in 
a PCEB may be a non-degenerate star or another compact 
object, e.g., binary pulsars (NS+NS), ultra-compact X-ray 
binaries (WD+NS) or AM CVn binaries (WD+WD). 

There is a vast literature on the theory of the common 
envelope phase a nd many published models of the popu- 
lation of PCEBs dPostnoy fc Yungelsonll2006l ; llben fc Liviol 
1 19931 ; iTaam fc Sandquistl l200ol v" Considerable debate con- 
tinues over fundamental aspects of the CE phase, much 
of it focused on the definition and value of the parameter 
qce = ^1°^ , where A_E or b is the change in orbital en- 
ergy of the binary during the CE phase and -E onv is the 
change in the binding energy of the red giant envelope. 
In theory, the outcome of the CE phase for a given bi- 
nary can be calculated given a known value for the param- 
eter qce- In practice, it is unclear whether the definition 
of E env should include energy sources othe r than the grav- 
itational potential energy of the envelope (ISoker fc Harpazl 
20031; iDewi fc Taurisll2000l; lHan et al.llioO^ ; iBeer et alj|2007l ; 



Webbinkll2008l ) . w hether evolution prior to Roc he lobe over- 
flow is important (|Eggletonl [20021 ; ISoked I2004T I or whether 
a bala nce of energies gives a valid description of the CE 
phase (|Nelemans et alJI^OOOlTNelemans fc Toutl2005l ). These 
debates are driven by the discovery and characterization 
of PCEBs whose properties, either individually or as a 
group, are not easily explained by a simple model of the 
CE ph ase, e.g., detached W D+WD "double d egenerate" bi- 
naries dMaxted et al.ll2002al lbl) and sdB stars (jMaxted et alj 
l2004l ). lDavis et al.l l 2009l ) have recently compared the results 
of binary population synthesis models for a range of input 
parameters and using different parametrizations of the com- 
mon envelope phase to the observed properties of 35 PCEBs. 
They find that the standard qce parametrization can ac- 
count for the observed properties of PCEBs with late-type 
companions to (pre-)white dwarfs, but cannot explain IK 
Peg which has an A8 companion star. They conclude that 
"the detection of more PCEBs with early type secondaries 



may shed further light on the CE phase" . They also find that 
there is a sharp decline in the observed number of PCEBs 
with periods > 1 day that cannot be reproduced by any 
of the models considered, even if the selection biases against 
longer period systems are accounted for. This suggests there 
are some important features of the common envelope phase 
that are not included in the current models. 

One way to better understand the common envelope 
(CE) phase is to compare the predictions of binary popu- 
lation synthesis models to the properties of a large, homo- 
geneously selected sample of PCEBs for which the selection 
effects are well understood. This has the advantage of giv- 
ing a result that is valid for the general population of that 
type of PCEB. Comparison of models to one or two systems 
may lead to spurious conclusions if one of the systems has 
had a peculiar history, e.g., V471 Tau appears to be a nor- 
mal PCEB containing a white dwarf and a K-dwarf, but its 
membership of the Hyades enables us to infer that it has a 
peculiar evolutionary history and is probably the result of 
the evolution of a triple star system (|0'Brien et al.ll200ll ). 

The Sloan Digital Sky Survey (SDSS) makes it pos- 
sible to construct a well-defined sample of PCEBs that 
can be identified from their optical colours. This is par- 
ticularly so for non-interacting binaries containing a WD 
with M-dwarf companions (WD+M). These binaries have 
distinctive colours that enable them to b e efficiently iden- 
tified from the SDSS phot o metry alone, (|Augusteiin et aCI 
2008; [Schreiber et al. l l2007l ; lYanmi r2009) and there is also 
spectroscopy of m a ny such systems from the same survey 
jHeller et al. 1 120081 ; iRebassa-Mansergas et all 120071 ). There 
are also many CVs with SDSS spectra that have been tar- 
geted for observ ation from their SD SS colours or observed 
serendipitously (|Szkodv et alj 120071 ). The selection effects 
for CVs are large and depend strongly on their accretion 
properties, whose dependence on the other properties of the 
binary are themselves poorly understood. For this reason, 
detached PCEBs are a better choice than CVs as a sample 
for testing models of CE evolution for WD binaries. 

Using SDSS photometry alone it is only possible to de- 
tect unresolved non-interacting binaries if the companion 
to the white dwarf has a spectral type later than about 
M0, the exa ct limit depending on t he temperature of the 
white dwarf (|Augusteiin et alj|2008l ). ISchreiber fc Gansickei 
(2003) calculated the temperature of the coolest WD that 
can be detected from a U-band excess as function of spectral 
type for a main-sequence companion. The limits are approx- 
imately 45,000K, 28,000K and 9,000K for spectral types K5, 
M0 and M6, respectively. The limits will be similar for the 
SDSS or any survey based on optical photometry because for 
cooler WDs the main-sequence companion is much brighter 
than a non-accreting WD at all optical wavelengths. 

To find non-accreting white dwarfs with F-, G- and Re- 
type companions or to find cool WDs with M-type com- 
panions requires a surv ey at UV or soft X-ray wavelengths. 
lLandsman et al] (|l996l ) identified two late-type stars with 
hidden hot white dwarfs (56 Per and HR 3643) from follow- 
up observations of late-type stars with a 1565A excess in the 
TD-1 ultraviolet sky survey. These are Sirius-like binaries, 
i.e., non-interacting systems with large orbital separations. 
Many more Sirius-like binaries were identified a s a result of 
the EUVE and ROSAT surveys (lBurleighll2002r). These sur- 
veys also identified several new PCEBs (|Burleigh et al.lll997l ; 
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Figure 2. Colour-colour diagrams for our sample of stars with reliable optical and ultraviolet photometry. The limit for selecting stars 
with an ultraviolet excess is shown in the left panel as a dashed line. Stars above this line are plotted with dots in the right panel. For 
clarity, we show the density of points in the right panel with contours rather than individual points. An approximate indication of spectral 
type for single stars based on their g — r colour is indicated in the bottom of the right panel. Stars for which we have obtained spectra are 
circled and stars with variable radial velocities are marked with diamonds. The dash-dot line with squares shows the combined colours 
of a white dwarf with T e fj =48,000K plus cool companions with spectral types from G2 to M5 (the blue-most point is the white dwarf 
alone). 



iBarstow et "all 1 1994 IVennes et all 1 19981 ). There have also 
been serendipitous discoveri es, e.g., from IUE spectroscopy 
|Boehm-Vitenselll993l . 119921 1. These serendipitous discover- 
ies can provide useful evidence for the existence of a partic- 
ular type of binary star, but it is difficult to know whether 
they should be included in any test of a binary population 
synthesis model. 

The advent of the GALEX all-sky UV survey 
l|Martin et alj|2005h makes it feasible to conduct a large sur- 
vey to detect hot white dwarfs in unresolved binaries with 
G-type and K-type companions and cool WD with early-M 
type companions. In this paper we report the first results of 
our survey to obtain follow-up spectroscopy of solar-type 
stars with a far-UV excess identified using GALEX and 
SDSS photometry. The main aim of the survey is to pro- 
duce a well-defined sample of detached PCEBs that cannot 
be identified using SDSS photometry alone. The properties 
of these PCEBs will be a strong test of population synthe- 
sis models. We also expect that this survey will find other 
interesting objects, e.g., CVs and other interacting binary 
stars. 



2 TARGET SELECTION 

To discover PCEBs in which the non-degenerate companion 
dominates the optical flux we require a sample of stars with 
reliable optical and ultraviolet photometry. We also require 
that the sources are bright so that the follow-up of the tar- 
gets can be done efficiently. To produce this sample of stars 
we used optical SDSS photometry and ultraviolet GALEX 
photometry. We started with the table provided by Multi- 
mission Archive at STScI of objects from SDSS DR4 cross 



matched within 4arcsec of objects in GALEX GR2 (table 
xSdssDr4) . From this table we selected objects satisfying the 
following criteria: 

• < 0.6° from the centre of the GALEX field of view; 

• the closest SDSS object to the GALEX source position; 

• GALEX magnitudes n < 25 and / < 25; 

• SDSS ^-magnitude g < 16.5; 

• SDSS morphological type class 6 (STAR); 

• "clean" SDSS photometry^] 

• proper motion > 0.01 arcsec/yr. 

The g magnitudes used in this paper are taken from the 
SDSS table entry psfMag_g, and similarly for other SDSS 
magnitudes. The resulting table has 2960 sources. Although 
the matching radius is 4arcsec, the positions for GR2 and 
DR4 agree to within 1 arcsec for 80% of the sources. The 
restriction on the source position in the field of view is to 
avoid spurious detect ions that are common in these regions 
(|Agiieros et al.ll2005T ). The proper motion information was 
taken from a table provided by SDSS DR5 and is based on 
the astrometry from the SDSS plus recalibrated USNO-B as- 
trometry. Objects detected at fewer than 6 epochs in SDSS 
plus USNO-B were excluded. This restriction on the proper 
motion was included to reduce the contamination of the sam- 
ple by background sources (galaxies and quasars). Without 
this criterion the sample size would be 8760 sources. The 
disadvantage of adding this proper motion criterion is that 
it introduces a kinematical bias into the survey. A complete 
treatment of this bias is beyond the scope of this paper. 



1 See http://skyserver.sdss.org/dr2/en/help/docs/ realquery. asp # flags 
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The majority of objects in this sample of 2960 sources 
are expected to be moderately bright stars with reliable op- 
tical and ultraviolet photometry. There are no stars brighter 
than g — 13 in the sample because stars brighter than this 
appear saturated in the SDSS images so they do not have 
reliable SDSS photometry. 

The far-UV (FUV) and near-UV (NUV) passbands of 
the GALEX instrument are shown in Fig. Q] and are com- 
pared to the UV spectra of three WDs and a very active 
GOV star. It can be seen that hot WDs (DA2 or hotter) 
can be easily distinguished from an active solar-type stars 
at FUV wavelengths. 

NUV and u-band fluxes can be strongly affected by 
chromospheric activity in late- type stars (Fig.[l|, so we used 
the f — g colour to identify stars with a UV excess. The dis- 
tribution of these sources in the / — g v. g — r and u — g v. 
g — r colour-colour diagrams are shown in Fig. [2] There is 
a well defined main-sequence relation in the optical colour- 
colour diagram, showing that the optical photometry for the 
majority of stars in our sample is reliable. Our criteria for 
selecting stars with a UV excess as shown in this figure are 

(/-(?)< 13( 5 -r) + 1.9 and (/-<?)< 4. 

There are 90 targets that satisfy these selection criteria. 



2.1 Cross-identifications for target objects 

We have used the SIMBAD database to search for cata- 
logued objects matching our target positions. The results 
are given in Table [1] and some individual objects are dis- 
cussed below. These stars give some impression of the type 
of object that may be discovered using our survey. In most 
cases, but not all (e.g. U Sex), we avoided observing objects 
where the nature of the star was already clear from existing 
observations. 

K UV 03134-0001 The s pectral type of this star given 
by IWegner fc Dupuisl (ll993T ) is "cont. + dM - Composite 
spectrum; shows blue continuum plus TiO features in the 
red". 

PQ Gem A well-kno wn intermediate polar l|Evans et al.l 
12009 ; iRosen et al.ll 19931 ). 

J0852+0313 The SDSS spectrum of this object shows 
that it is a QSO with a redshift z = 0.297. There is a faint 
companion within a few arcseconds to this QSO visible in 
the SDSS images that is unresolved in the DSS images from 
which the USNO-B positions are measured. We conclude 
that the proper motion value for this objec t is spurious. 

BZ UMa This is a dwarf nova (CV, iNeustroev et all 
120061 ) . 

PG 1056+324 I Wagner et"all (|l988l ) give a spectral type 
of "sdB" for this star. 

BE UMa This is a well-studied detached PCEB con- 
taining a very hot (pre-)white dwarf and a low-mass, 
non-degenerate comp anion with an orbital period of 2.3 d 
^Ferguson et alJll999h . 

SDSS J14091 6. 11+382832.1 The a utomatic classifica- 
tion algorithm of lEisenstein et all l|2006T ) identifies this star 
as a hot subdwarf. Visual inspection of the SDSS spectrum 
shows hints of absorption lines due to a late-type companion 
(G-band, Mg-b, Call IR triplet). 



SBS S 1422+497 The spectral type given by IStepanianl 
l|2005l ) is "DA". 

PG1657+416 This is a pul sating sdB star w ith a G5 
main-sequence companion star |Oreiro et alj|2007h . 

2.2 Observations 

We obtained spectroscopy for 36 targets from our sample us- 
ing the Intermediate Dispersion Spectrograph (IDS) on the 
2.5m Isaac Newton Telescope. We used the H2400B grat- 
ing with a 1.2arcsec slit and the EEV10 charge coupled 
device (CCD) detector to obtain spectra with a resolution 
of 0.45A sampled at 0.22A/pixel. The unvignetted portion 
of the CCD covers the wavelength range 4875 -5375 A. We 
also extracted the region of the spectrum around the H/3 
line that lies in the vignetted region of the CCD. Spec- 
tr a were extrac ted using the optimal extraction algorithm 
of iHorn el jl986). Observations of each star were bracketed 
with arc spectra and the wavelength calibration established 
from these arcs interpolated to the time of mid-exposure. 
We obtained a single spectrum of 19 of our targets on 2007 
Jan 29. The majority of the spectra were obtained in the 
interval 2008 Dec 20 to 2008 Dec 26. The average spectrum 
of each target is shown in Fig. [3] 

We measured the radial velocity (RV) of our targets 
using cross correlation against a spectrum of t he RV stan- 
dard star HD 65583 (G8V, V r = +14.85 kms" 1 ; lUdrv et ail 
1999). The spectral range used for the cross-correlation was 
4868-5395A. We obtained spectra of the RV standard star 
HD 3765 on every night of the 2007 Jan and 2008 Dec runs 
to monitor the stability of our RV measurements. The stan- 
dard deviation of the RV measurements from these spectra 
is 2.79 kms -1 . This is much larger than the precision with 
which the peak of the cross-correlation function (CCF) can 
be located and is likely to be dominated by the motion of 
the stellar image within the spectrograph slit. We added this 
estimate of the external error in quadrature to the internal 
error of the RV measurements of our targets to obtain a to- 
tal standard error estimate for each RV measurement. The 
exposure times used for our targets (300 - 1800s) are much 
longer than those used for HD3765 (30 -60s). Guiding er- 
rors during a longer exposure will smear the light of the star 
across the slit so the external error due to image motion in 
the slit is likely to be less for our targets than for HD 3765, 
i.e., the total standard errors we have adopted are likely to 
be pessimistic. 

For the set of RV measurements of each target V r ,i; i = 
1, . . . , JVrv, we calculate the weighted mean RV value, V r . 
To test for variability of the RV we calculate the probability, 
p, of obtaining the observed value of \ 2 or greater for the 
model V r ,i = V r assuming normally distributed errors, i.e., 
a small value of p indicates that the RV measurements are 
not constant. The value of p for each target are also given 
in Table [2] as log p. 

Also shown in Table [2] is a crude estimate of the equiva- 
lent width of the H/3 line, EW(H/3) in the average (median) 
spectrum of the star after shifting the individual spectra to a 
common wavelength scale according to their measured RVs. 
The equivalent width was calculated by numerical integra- 
tion in a window ±5A wide around the rest wavelength of 
the line. Negative values of EW(H/3) indicate that the H/3 
line is in emission. 
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Figure 3. Spectra of our targets ordered by their g — r colour. Spectra have been co-added on a common wavelength scale, normalized, 
smoothed and re-binned for display and are offset by 1 for clarity. Spectra of some radial velocity standard stars with known spectral 
types as shown in the label are included for comparison. The rest wavelength of the H/3 line is indicated with a dotted line. 



3 RESULTS 

3.1 Detection efficiency 

For each star we calculate the probability of detecting a com- 
p anion from the r adial velocity variations using the method 
of iMaxted etTail i|200ll) . We assume a mass of 0.75Mq for 



the visible star because this is likely to be close to the av- 
erage mass of our target stars. We assume a mass ratio of 
q — 0.75 because this value then gives a white dwarf com- 
panion mass of 0.56-Mq, which is a typical mass for a white 
dwarf. The results for two stars are shown in Fig. [4] We find 
that the average detection efficiency near a period of 10 d is 
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Table 1. Cross-identifications for objects from our sample with sources listed by SIMBAD. 



ID 




^2000 


<52000 




9 


u- g 


g-r 


f-9 


Other ID 


Notes 


J0316- 


-0009 


03 16 


00.2 


+00 


09 


47 


16.45 


0.52 


0.65 


0.06 


KUV 03134-0001 




J0751- 


L 1444 


07 51 


17.3 


+14 


44 


24 


14.18 


—0.06 


0.18 


1.44 


PQ Gem 


Intermediate polar 


J0824- 


r3028 


08 24 


34.1 


+30 


28 


55 


15.25 


0.09 


—0.07 


—0.44 


PG 0821+306 


QSO 


J0852- 


r0313 


08 52 


59.2 


+03 


13 


21 


16.22 


0.08 


0.06 


1.30 


2MASS J08525922+0313207 




J0853- 


-5748 


08 53 


44.2 


+57 


48 


41 


16.40 


—0.41 


0.35 


0.74 


BZ UMa 


Dwarf nova 


J0941- 


-4706 


09 41 


55.2 


+47 


06 


38 


16.12 


2.06 


1.28 


3.89 


XMS J094154. 6+470637 




J0957- 


r0340 


09 57 


25.4 


+03 


40 


05 


14.91 


1.18 


0.25 


2.61 


U Sex 


RR Lyrae star 


J1014- 


r0033 


10 14 


54.9 


+00 


33 


37 


16.27 


0.08 


0.16 


0.66 


QSO B1012+008 




J1042- 


-6442 


10 42 


34.8 


+64 


42 


05 


15.00 


1.57 


1.28 


3.67 


1RXS J104237. 0+644217 




J1054- 


L 3600 


10 54 


50.6 


+36 


00 


13 


16.48 


2.47 


0.96 


3.12 


LEDA 3088408 




J1059- 


r 3206 


10 59 


05.2 


+32 


06 


21 


14.64 


0.14 


—0.09 


—0.45 


PG 1056+324 


sdB 


J1121- 


r0012 


11 21 


19.7 


+00 


12 


12 


14.90 


1.18 


0.42 


3.38 


2MASS J11211965+0012120 




J1155- 


-0640 


11 55 


28.1 


+06 


40 


18 


15.70 


1.05 


0.43 


1.91 


Pul-3 850246 




J1157- 


r 4856 


11 57 


44.8 


+48 


56 


18 


15.78 


-0.28 


-0.10 


-1.72 


BE UMa 


PCEB, P=2.3d 


J1300- 


h4139 


13 00 


40.1 


+41 


39 


37 


15.07 


0.94 


0.34 


3.05 


BPS BS 16938-0017 




J1409- 


h3828 


14 09 


16.1 


+38 


28 


32 


16.03 


-0.13 


-0.16 


-1.03 


SDSS J140916. 11+382832.1 




J1418- 


r3705 


14 18 


31.8 


+37 


05 


41 


15.41 


1.27 


0.46 


2.99 


2MASS J14183182+3705411 




J1424- 


^4929 


14 24 


40.5 


+49 


29 


58 


16.24 


0.10 


0.49 


-1.06 


SBSS 1422+497 


DA 
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41.8 
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16 
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PG 1657+416 


sdB+G5V 




0.0 L 
0. 



log(P/d) 



100.0 



Figure 4. The fraction 
assuming a mass ratio q - 



of binaries detected with logp < —3 
= 0.75 and a mass for the visible star of 
0.75Mq using data with the same sampling and standard errors 
as that obtained for J0801+2143 (solid line) and J0909+4521 
(dashed line). The detection efficiency in each bin is the average 
of 100 periods evenly distributed between the bin limits. 



79% and that the detection efficiency is more than 50% for 
periods of less than 10 d for all 36 of the stars observed. 

3.2 Notes on individual objects 

J0754+2035 We obtained a single spectrum of this star 
with the IDS and the H2400B grating at a central wave- 
length of 4200A. This spectrum shows strong emission in 
the cores of the Call H and K lines (EW » -3.5A) for 
this star (Fig. [5}. The width of the CCF in this star shows 
this star is a rapid rotator and so is expected to be chromo- 
spherically active. There is some weak evidence for a small 
amount of RV variability in this star (« lOkms -1 ), which 
may be a result of star spots rather than binarity. 

J0800+5231 There is an offset of about 10km s" 1 be- 
tween the RVs measured for this star in 2007 and 2008, but 




J0800+5231 



J0735+2035 




3950 



4000 



4050 



4100 



Wavelength (A) 



Figure 5. Spectra of selected targets in the region of the Call 
H and K lines. The spectra have been normalized, smoothed and 
vertically offset for clarity. 



no significant variability over 4 nights from the 2008 data 
alone. This may be a binary with a long orbital period. A 
single spectrum near 4200A obtained with the IDS (Fig. [5]) 
suggests that the spectral type is late-G, which is consistent 
with the u — g and g — r colours of the star. 

J0814+5328 The spectra near the H/3 line are noisy but 
the absorption line does appear to be shallow compared to 
other stars of the same colour/spectral-type (Fig. |3). The 
absorption line may be "filled-in" by weak emission at this 
wavelength. 

J0851+2903 This is a similar case to J0754+2035. A sin- 
gle spectrum around 4200A shows Call H and K lines in 
emission (EW « — 3A, Fig. O, and there is some evidence 
of RV variability at the level of 20km s _1 which may be the 
result of star spots in this rapidly rotating star. 

J0905+3256 This star has a large heliocentric RV (« 
— 320kms~ 1 ) so is likely to be a halo star. 

J0941+47 06 This star is an X- ray source detected by the 
XMS survey (jBarcons et al.|[2007h . The RV changes by more 
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Table 2. Summary of our radial velocity measurements. N rv is the number of radial velocity measurements. At is the baseline 
of the observations. See text for the definition of logp. FWHM is the mean full-width at half-maximum of the CCF. Where a * 
symbol appears in the final column further remarks on the star appear in Section 13.21 



ID 


RA 

(J2000) 


Dec 
(J2000) 


9 


u - 


-9 


g-r 


f-9 




logp 


At 

(days) 


EW(H/3) 

(A) 


FWHM 

(kms" 1 ) 




J0357— 0655 


03:57:16.0 


—06:55:04 


16.07 


1 


.75 


1.00 


1.58 


5 


—0.10 


693.1 


1.4 


70 




J0735+2035 


07:35:34.5 


+20:35:32 


15.87 


2 


.04 


0.84 


3.91 


5 


—0.51 


696.1 


0.7 


67 




J0740+2319 


07:40:33.6 


+23:19:26 


14.56 


1 


.53 


0.59 


2.21 


2 


—0.02 


691.1 


1.8 


65 




J0751+1940 


07:51:32.0 


+19:40:55 


16.46 


1 


.72 


0.88 


1.66 


4 


0.00 


691.2 


1.5 


69 




J0754+2035 


07:54:07.5 


+20:35:52 


16.12 


2 


.13 


1.08 


3.36 


6 


-1.12 


696.1 


-0.4 


113 


* 


J0754+2847 


07:54:04.2 


+28:47:01 


16.10 


1 


.67 


0.62 


3.99 


2 


—0.25 


691.2 


1.4 


67 




J0755+2424 


07:55:41.9 


+24:24:36 


16.27 


1 


.14 


0.35 


3.55 


2 


-0.35 


692.1 


2.4 


86 




J0800+5231 


08:00:38.6 


+52:31:29 


16.28 


1 


.30 


0.56 


3.20 


i 


-2.69 


696.0 


1.1 


69 


* 


J0801+2143 


08:01:07.1 


+21:43:37 


15.09 


1 


.85 


0.74 


3.70 


2 


-0.17 


691.2 


1.2 


69 




J0807+2225 


08:07:36.5 


+22:25:37 


16.21 


2 


.02 


0.85 


2.72 


2 


—0.03 


691.2 


1.0 


73 




J0814+5328 


08:14:43.4 


+53:28:05 


16.48 


1 


. 13 


0.62 


2.28 


2 


-0.23 


692.0 


0.1 


92 


* 


J0816+3707 


08:16:07.1 


+37:07:04 


16.39 


1 


.67 


1.23 


2.57 


2 


-0.11 


691.1 


0.9 


69 




J0824+2512 


08:24:56.5 


+25:12:47 


15.83 


1 


.65 


0.64 


3.36 


2 


—0.25 


691.2 


1.7 


67 




J0834+4104 


08:34:50.0 


+41:04:23 


15.73 


1 


.39 


0.59 


1.24 


2 


—0.14 


691.2 


1.5 


65 




J0851+2903 


08:51:37.2 


+29:03:30 


16.24 


1 


.65 


1.05 


1.66 


9 


—2.16 


696.1 


-1.3 


120 




J0858+0512 


08:58:03.2 


+05:12:49 


15.84 


2 


.11 


0.99 


3.03 


2 


-0.11 


691.3 


1.1 


69 




J0902+3847 


09:02:04.1 


+38:47:38 


14.74 





.90 


0.24 


1.51 


6 


-0.99 


693.2 


2.1 


76 




J0903+0949 


09:03:51.2 


+09:49:35 


16.38 





.98 


0.25 


3.74 


3 


-0.09 


4.1 


2.7 


121 




J0905+3256 


09:05:24.4 


+32:56:02 


15.24 





.96 


0.27 


3.83 


1 


—0.37 


4.0 


2.1 


76 


* 


J0909+4521 


09:09:38.0 


+45:21:42 


15.63 


2 


.08 


0.90 


3.37 


3 


—0.21 


4.1 


0.5 


71 




J0915+5914 


09:15:42.7 


+59:14:57 


15.84 


2 


.01 


1.05 


2.69 


3 


-0.09 


4.0 


0.7 


74 




J0917+4442 


09:17:10.3 


+44:42:26 


15.84 


1 


.95 


0.99 


2.50 


1 


-0.07 


4.0 


0.6 


69 




J0923+5735 


09:23:23.1 


+57:35:40 


15.25 


1 


.05 


0.48 


1.22 


3 


—0.12 


4.1 


0.4 


60 




J0926+2811 


09:26:16.6 


+28:11:13 


15.98 





.92 


0.36 


3.90 


4 


—0.03 


4.1 


1.6 


61 




J0941+4706 


09:41:54.6 


+47:06:38 


16.12 


2 


.06 


1.28 


3.89 


4 


< -40 


1.0 


-1.5 


96 


* 


J0957+0340 


09:57:25.4 


+03:40:05 


14.91 


1 


.18 


0.25 


2.61 


3 


—0.36 


1.1 


4.3 


76 




J0959+3621 


09:59:02.4 


+36:21:03 


14.49 


1 


.36 


0.50 


3.39 


2 


-0.07 


1.0 


2.3 


58 




J1016+4935 


10:16:57.1 


+49:35:14 


16.11 


1 


.55 


1.19 


1.44 


i 


< -40 


3.9 


-0.8 


71 


* 


J1019+3821 


10:19:37.4 


+38:21:12 


15.41 


2 


.19 


0.96 


3.05 


2 


-0.46 


1.0 


0.9 


68 




J1037+4200 


10:37:08.3 


+42:00:39 


15.78 


2 


.11 


1.10 


2.69 


2 


-0.07 


1.0 


0.3 


69 




J1042+6442 


10:42:34.8 


+64:42:05 


15.00 


1 


.57 


1.28 


3.67 


9 


< -40 


3.0 


-5.1 


180 


* 


J 1059+3206 


10:59:05.2 


+32:06:21 


14.64 





.14 


-0.09 


-0.45 


2 


-1.06 


2.9 


2.5 


66 




J2104+0026 


21:04:01.4 


+00:26:54 


16.29 


1 


.53 


0.56 


3.38 


2 


-0.05 


4.0 


2.7 


57 




J2113+1107 


21:13:16.2 


+11:07:48 


16.31 


1 


.17 


0.53 


1.85 


2 


-0.24 


4.0 


1.5 


70 




J2320-1045 


23:20:40.4 


-10:45:11 


15.22 


1 


.32 


0.50 


3.30 


6 


-0.02 


5.0 


1.6 


63 




J2347+1335 


23:47:42.4 


+13:35:57 


15.54 


2 


.38 


1.07 


3.89 


1 


-0.01 


5.0 


0.6 


69 





than 50km s _1 in less than 3 hours. This is too large to be 
explained by star spots. The H/3 emission line has a similar 
width to the CCF and shows the same RV variation, i.e., the 
H/3 emission is due to irradiation of and/or chromospheric 
activity on the M-dwarf and is not due to an accretion disc. 
This is almost certainly a PCEB. 

J0957+03 40 This is an RR L yr-type variable star of 
type RRab (jKukarkin et all ll97lT ). The amplitude of the 
variability of these st ars at UV wavelengths is very large 
|Wheatlev et al.ll2005l ). so it is not surprising that optical 
and UV photometry obtained at different epochs result in a 
peculiar position for the star in the UV-optical colour-colour 
diagram. 

J1016+4935 The RV of this star changes by » lOOkms" 1 
over 3d. The H/3 emission line has a similar width to the 
CCF and shows the same RV variation, i.e., the H/3 emission 
is due to irradiation of and/or chromospheric activity on the 
M-dwarf. This is almost certainly a PCEB. We have too few 
data to estimate the orbital period. 



J1042+6442 This star is an X-ray source de tected by 
the ROSAT-FSC survey (jMickaelian et alJ liooril. The RV 
changes by approximately 250km s -1 in less than 3 hours. 
Additional spectroscopy has been obtained and is discussed 
in the next section. 



3.3 J1042+6442 

The data for this star obtained during 2008 Dec clearly show 
that the RVs vary with a period ~ 0.2 d. We obtained fur- 
ther spectroscopy with the same instrument and configura- 
tion during Jan 2009 to measure the orbital period and mass 
function of this star. A further 17 spectra were obtained, al- 
though some of these have poor signal-to-noise due to poor 
seeing. The RV of J 1042+6442 was measured as above for all 
the spectra and are given in Table [3] We then searched for 
periodicity in these RV measurements by fitting sine waves 
at a range of trial frequencies. The period is found to be 
close to 0.19767 d and is unambiguous (Fig. [6}. Also given 
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Figure 6. Periodogram of our RV measurements for J1042+6442. 
The value of x 2 is shown for a least-squares fit of a sine wave at 
each frequency. 
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Figure 7. Upper panel: RV measurements of J1042+6442 with 
the circular orbit fit by least-squares plotted as a function of or- 
bital phase for P = 0.19767 d. Lower panel: Equivalent width 
of the H/3 emission line. . Measurements with large uncertainties 
have been excluded for clarity. 

in Table [3] is the equivalent width of the H/3 emission line 
measured by numerical integration in a window ±250kms _1 
around the position of the line based on the measured RV. 
We used a least-squares fit of a circular orbit to the RV 
measurements to derive the parameters shown in Table 2] 
The scatter of the RV measurements around the best fit is 
larger than expected given the standard errors of the RV 
measurements. We account for this by adding a quantity 
cr cx in quadrature to the standard errors of the RV mea- 
surements prior to calculating the fit. The physical origin 
of ow may be related to the extreme chromospheric activ- 
ity that is expected for such a rapidly rotating M-dwarf. 
The RV measurements, circular orbit fit and EW(H/3) mea- 
surements are shown in Fig. [7] The H/3 emission line has a 
similar width to the CCF and shows the same RV variation 
so it originates from the M-dwarf. There is a clear cosine- 
like variation of EW(H/3) with orbital phase, as expected if 
it is due to the irradiation of one side of the M-dwarf by a 
hot companion. We applied a rotational broadening function 



to the spectra of HD 79210 (MOV) for various values of the 
projected equatorial rotational velocity, V ro t sini where i is 
the inclination, and measured the FWHM of the resulting 
CCF. From this calibration we estimate that the FWHM of 
180km s" 1 observed for the CCF of J1042+6442 corresponds 
to V r ot sin i~ 90km s _1 . This corresponds to a "projected 
radius" Rsini ~ 0.36-Rq. We have compa red the colour of 
this st ar to the colour-colour diagrams of lAugusteiin et al.l 
(2008). From the (r — i) v. (i — z) diagram (their Fig. 2) we 
estimate a spectral type of ~M1 for J1042+6442. This spec- 
tral type is consistent with the appearance of the spectrum 
shown in Fig. [3] e.g., the weak TiO bandhead near 4950A 
is slightly stronger than the MOV star HD 79210. The mass 
of an M1V star is about 0.5Mq and the radius is about 
0.5-Rq. Combined with the mass function f m — 0.093Mq, 
we estimate a system inclination i ~ 45° so the mass for the 
unseen companion to the Ml star « 0.75Mq. There is large 
scatter in the relationship between mass, radius and spectral 
type for M-dwarfs, but these estimates are quite consistent 
with the companion being a white dwarf. 

The radius of the Roche lobe in this system is ~ 
0.5-Rq, which is comparable to the radius expected for 
an Ml-type dwarf star. This raises the possibility that 
there is accretion onto the companion from the M-dwarf 
through the inner Lagrangian point. However, there is no 
evidence for an accretion disc or other accretion struc- 
tures from the H/3 emission line. J1042 +6442 has b een de- 
tected in the ROSAT All Sky Survey (|Vogesll200ol ) with a 
count rate of 0.046 ± 0.011 cts/s and with a hardness ra- 
tio HR1 = —0.29 ± 0.20, implying a relatively soft source. 
The galactic column density in the direction of J1042+6442 
is 1.2 x 10 20 cm -2 (jDickev fc Lockmanl Il99(ih . Given the 
short distance, 67 pc, we assume that about half of that 
column is in front of J1042+6442. Adopting a thermal 
Bremsstrahlung spectrum and a plasma temperature in 
the range 1 — lOkeV results in an X-ray luminosity of 
J1042+6442 of L x = (2 - 6) x 10 29 erg/s, which is roughly 
an order of magnitude lower t han the typical X-ra y lumi- 
nosity of cataclysmic variables l|Verbunt et al.lll997t ). but is 
consistent with the X-ray luminos ity of rapidly rotating low- 
mass stars l|Pizzolato et al.ll2003l ). More specifically, adopt- 
ing a 0.5 keV Raymond-Smith plasma, appropriate for coro- 
nal emission, results in Lx = 2 x 10 29 erg/s, which is very 
similar to the X-ray luminosities of the eclipsing short-period 
low-mass stars GUBoo, YYG em, and CU Cnc (lRibasll2003l ; 
lLopez-Morales fc Ri bas 2005). We conclude that the X-ray 
emission from J1042+6442 is consistent with coronal emis- 
sion from the main-sequence component, corroborating the 
detached pre-CV nature of the system. 

3.4 White dwarf effective temperatures 

We have estimated the effective temperature of the white 
dwarfs in the three new PCEBs we have identified based 
on the available SDSS and GALEX photometry. We first 
use th e SDSS i — z colou r of the system and the calibra- 
tion of ICovev et alj l|2007h to estimate the spectral type of 
the dwarf star in these binaries. We then use the observed 
z-band magnit ude togethe r with the Mj and z — J cali- 
bration from I Covey et al.l to estimate the distance to the 
dwarf star. J-band magnit udes are taken from the 2MASS 
catalogue (|Skrutskiel l2006V The contribution of the white 



© 2008 RAS, MNRAS OOO.mfiOl 



A survey for PCEBs using GALEX and SDSS 9 



Table 3. Radial velocity and EW(H/3) measurements for 
J1042+6442 



MJD 



V r 
(kms" 1 ) 



-EW(H/3) 
(A) 



54823 


2846 


112.0 ± 3 


7 


4.38 ± 0.27 


54826 


1419 


-163.7 ± 2 


9 


3.78 ± 0.24 


54826 


2602 


91.4 ± 1 


8 


5.48 ± 0.24 


54826 


2638 


78.3 ± 3 





5.52 ± 0.25 


54826 


2887 


-64.9 ± 2 


6 


5.29 ± 0.26 


54826 


2923 


-93.3 ± 2 


8 


4.80 ± 0.25 


54826 


2959 


-119.4 ± 2 


i 


5.59 ± 0.27 


54826 


2995 


-91.5 ± 2 


9 


5.50 ± 0.27 


54826 


3032 


-117.0 ± 2 


6 


5.55 ± 0.29 


54843 


1881 


26.2 ± 5 


3 


5.87 ± 0.67 


54843 


2356 


142.9 ± 8 


2 


7.48 ± 1.0 


54843 


2643 


45.8 ± 15 







54843 


2796 


-57.9 ± 8 


6 


5.32 ± 1.1 


54843 


9713 


15.6 ± 3 


6 


2.09 ± 0.36 


54844 


0278 


184.3 ± 4 


7 


3.50 ± 0.21 


54844 


1461 


-76.2 ±3 


7 


3.20 ± 0.25 


A QA A 


zolo 


(0.0 ± 


U 


A QA I r. r o 
4.04 ± U.OO 


54844 


3147 


-147.1 ±6 


6 


1.16 ± 0.72 


54844 


9890 


155.1 ± 2 


i 


2.04 ± 0.28 


54845 


0185 


145.3 ± 2 


5 


4.01 ± 0.21 


54845 


0760 


-79.2 ± 2 


8 


4.84 ± 0.32 


54845 


1065 


-166.5 ±4 


8 


3.26 ± 0.20 


54845 


2039 


172.0 ± 2 


9 


3.39 ± 0.20 


54845 


2304 


118.9 ± 2 


7 


3.76 ± 0.20 


54845 


2375 


77.6 ± 2 


3 


4.98 ± 0.22 


54845 


2774 


-134.4 ±4 


1 


3.17 ± 0.37 



Table 4. Circular orbit least-squares fit to our RV measurements 
of J1042+6442. The model is V r = 7 + Xsin[27r(T-T )/P] where 
T is the time of mid-exposure. The effect of orbital smearing due 
to the finite exposure time have been included. The quantity <r e x 
is added in quadrature to the standard errors in Table [3] prior to 



calculating the least-squares fit and is chosen so that 
the number of degrees of freedom in the fit. 



Parameter 



Value 



Period (days) 
To (HJD) 
7 (kms -1 ) 
K (kms" 1 ) 
o- ex (kms -1 ) 

JVdf 
,^2 



0.197669 ± 0.000022 
2454834.7820 ± 0.0012 
-0.1 ±4.4 
166.7 ± 5.3 
18 
22 
21.8 



dwarf at these wavelengths is negligible. Synthetic spectra 
fo r the white dwarf ge nerated using the models described 
in iKoester et all (2005) were then matched "by-eye" to the 
observed GALEX FUV and NUV fluxes by varying the ef- 
fective temperature assuming the white dwarf has a sur- 
face gravity logy = 8.0. The results are given in Table [S] 
It is difficult to make an accurate estimate of the uncer- 
tainty in these estimates. For example, the three FUV fluxes 
reported for J1042+6442 in GALEX GR4 are 139 ± 11, 
112.2 ± 2.7 and 133.1 ± 3.9, which are clearly not consis- 
tent with a single value. We also found differences of 10- 



Table 5. Estimate of the white dwarf effective temperature 
(T G ff wd) based on the SDSS and GALEX for our three newly 
discovered PCEBs. The spectral type (SpTy) and distance (d) of 
the dwarf star based on the observed i — z, colours and z magni- 
tudes are given in columns 2 and 3, respectively. 



Star 


SpTy 


d 


T e ff ,WD 


i 


i — z 


z — J 






(pc) 


(K) 








J0941+4706 


K7V 


280 


12000 


14.05 


0.31 


1.19 


J1016+4935 


M2V 


190 


14400 


13.96 


0.36 


1.22 


J1042+6442 


M3V 


67 


9800 


12.26 


0.54 


1.33 



20 percent when we c ompared the distances estimated us- 
ing the calibration of Covey et al. t o thos e derived using 
the calibrations of iDavenport et alj l|2006j ). Chan ging the 
assumed value of logy by ±0.5 changes the estimate of the 
effective temperature by about ±1000K for J0941±4706 and 
J1042±6442, and ±2000K for J1016±4935. Given these fac- 
tors, our opinion is that a fair estimate of the errors in the 
effective temperatures of the white dwarfs is a few thousand 
Kelvin. 



4 DISCUSSION 

The survey strategy we have adopted has identified one new 
PCEB (J1042±6442) and two stars that are also likely to be 
PCEBs (J0941±4706, J1016±4935) based on the amplitude 
and timescale of the RV variability and the narrow width of 
their H/3 emission lines. These new identifications are in ad- 
dition to two CVs (PQ Gem and BZ TJMa) and one pre-CV 
(BE UMa) in the sample that satisfy our selection criteria. 
Further observations will be required to determine the ba- 
sic parameters of J0941±4706 and J1016±4935. We have 
compared t he colours of these two st ars to the colour-colour 
diagrams of lAugusteiin et alj (|2008l ). The dwarf star domi- 
nates the flux at these wavelengths so the colours are indis- 
tinguishable from those of single stars. The other colours of 
J0941±4706 are also indistinguishab le from those of a sin- 
gle star according to the criteria of lAugusteiin et al.l . i.e., 
this is a PCEB that cannot be identified from its SDSS 
colours alone. The colours of J1016 ±4935 and J1042± 6442 
do fall within the selection criteria of lAugusteiin et"ai] in the 
(u — g) v. (g — r) plane (although they are near the limit), 
but not the (g — r) v. (r — i) plane. None of these stars 
appear in the catalogue of lAugusteiin et all because they 
exceed the brightness limit of g = 16.5 set by these authors 
to avoid stars with unreliable photometry. The effective tem- 
peratures we have estimated for the white dwarfs in these 
binaries (Section 13. 4|l are much lower than the temperature 
of the coolest WD that can be dete c ted fr om a U-band ex- 
cess given bv lSchreiber fc Gansickel l|2003l ). which supports 
our conclusion that the combination of GALEX and SDSS 
photometry can be used to identify PCEBs that are difficult 
or impossible to identify using SDSS photometry alone. 

Rather than setting a brightness limit, we have used 
the flags provided with the SDSS photometry to identify 
stars with reliable photometry. The brightness at which sat- 
uration of the CCD will produce unreliable photometry will 
depend on the filter under consideration, the colour of the 
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star and the seeing during the observation. It is noticeable 
that the SDSS photometry for several of our targets were 
obtained in worse-than-average seeing conditions. A full ex- 
ploration of how these factors affect our survey is beyond 
the scope of this paper. 

Two of the stars we have observed (J0754+2035 and 
J0851+2903) are very rapidly rotating K-dwarfs with strong 
chromospheric emission that show small shifts in their mea- 
sured RVs. These may be the result of wind-accretion in- 
duce d rapid rotation (WIRRing stars; Ijeffries fe Stevens! 
1996). The FUV-excess in these stars would then be due 
to the white dwarf remnant of the asymptotic giant branch 
(AGB) star from whose wind the K-dwarf accreted both ma- 
terial and angular momentum. The current separation of the 
WD and K-dwarf can be large in this scenario (~ lOOau) so 
it may be possible to directly observe these WD using high 
resolution imaging. The apparent RV variability of these 
stars is likely to be dominated by spurious shifts associ- 
ated with their extreme magnetic activity, i.e., star-spots, 
and is unlikely to be due to orbital motion. J0851+2903 can 
be identifi ed as a binary st ar containing a WD using the 
criteria of lAugusteiin et alj in the (g — r) v. (r — i) plane, 
but J0754+2035 would be indistinguishable from a single 
K-dwarf using the SDSS photometry alone. 

The majority of the stars in our sample appear to be 
normal late-type dwarfs. There is little contamination from 
unidentified quasars as a result of our proper motion se- 
lection. The FUV excess in the targets that do not show 
radial velocity shifts in our data may be due to a hot sub- 
dwarf of hot white dwarf companion to the late-type star 
in wide orbit. High resolution imaging or RV measurements 
with greater precision than those use here could be used to 
confirm this hypothesis. 



5 CONCLUSIONS 

We have shown that complementing SDSS photometry with 
GALEX photometry will enable us to identify new PCEBs 
and WIRRing stars that cannot be identified from SDSS 
photometry alone. The new binary stars we have identified 
with our follow-up spectroscopy contain late-K or early-M 
dwarf stars. A more complete survey of targets identified 
using the methods we have developed here will be required 
to show whether or not the lack of PCEBs containing dwarf 
stars earlier than late-K is a real feature of this population. 
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